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in an effective field theory description of nuclear forces. By using a toy 

model of the separable NN potential leading to logarithmic singularities in 
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which is governed by a generalized dynamical equation with a nonlocal-in-time 
interaction operator. 

PACS: 21.45.+V, 02.70.-c, 24.10.-i 

Keywords: Effective field theory; Nucleon dynamics; Nonlocal-in-time 
interaction. 



Understanding how nuclear forces emerge from the fundamental theory 
of quantum chromodynamics (QCD) is one of the most important problem of 
quantum physics. To study hadron dynamics at scales where QCD is strongly 
coupled, it is useful to employ effective field theories (EFT's) [1] being an 
invaluable tool for computing physical quantities in the theories with disparate 
energy scales. In order to describe low energy processes involving nucleons 
and pions, all possible interaction operators consistent with the symmetries 
of QCD are included in an effective Lagrangian of an EFT. However such a 
Lagrangian leads to ultraviolet (UV) divergences that must be regulated and a 
renormalization scheme defined. A fundamental difficulty in an EFT description 
of nuclear forces is that they are nonperturbative, so that an infinite series 
of Feynman diagrams must be summed. Summing the relevant diagrams is 
equivalent to solving a Schrodinger equation. However, an EFT yields graphs 
which are divergent, and gives rise to a singular Schrodinger potential. For 
this reason N-nucleon potentials are regulated and renormalized couplings are 
defined [2]. Nevertheless, a renormalization procedure does not lead to the 
potentials satisfying the requirements of ordinary quantum mechanics, and 
consequently after renormalization nucleon dynamics is not governed by the 
Schrodinger equation. This raises the question of what kind of equation governs 
nucleon dynamics at low energies. The Schrodinger equation is local in time, 
and the interaction Hamiltonian describes an instantaneous interaction. This is 
the main cause of infinities in the Hamiltonian formalism. In Ref. [3] it has been 
shown that the use of the Feynman approach to quantum theory in combination 
with the canonical approach allows one to extend quantum dynamics to describe 
the evolution of a system whose dynamics is generated by nonlocal in time 
interaction. A generalized quantum dynamics (GQD) developed in this way 
has been shown to open new possibilities to resolve the problem of the UV 
divergences in quantum field theory [3]. An equation of motion has been 
derived as the most general dynamical equation consistent with the current 
concepts of quantum theory. Being equivalent to the Schrodinger equation in 
the particular case where interaction is instantaneous, this equation permits the 
generalization to the case where the interaction operator is nonlocal in time. 
Note that there is one-to-one correspondence between nonlocality of interaction 
and the UV behavior of the matrix elements of the evolution operator as a 
function of momenta: The interaction operator can be nonlocal in time only in 
the case where this behavior is "bad", i.e. in a local theory it leads to the UV 
divergences. For this reason one can expect the nucleon dynamics that follows 
from renormalization of an EFT to be governed by the generalized dynamical 
eauation with nonlocal-in-time interaction ODerator. In the Dresent DaDer this 



problem is investigated by using the example of a toy model of the separable 
NN potential which gives rise to logarithmic singularities in the Born series. We 
show that the T matrix obtained in Ref. [4] by dimensional regularization of this 
model does not satisfy the Lippmann-Schwinger (LS) equation but satisfies the 
generalized dynamical equation with a nonlocal-in-time interaction operator. 
In the GQD the following assumptions are used as basic postulates: 
(i) The physical state of a system is represented by a vector (properly by 
a ray) of a Hilbert space. 

(ii) An observable A is represented by a Hermitian hypermaximal operator 
a. The eigenvalues a r of a give the possible values of A. An eigenvector \ip\f' > 
corresponding to the eigenvalue a r represents a state in which A has the value 
a r . If the system is in the state \ip >, the probability P r of finding the value a r 
for A, when a measurement is performed, is given by 

P r =<^|Pv;|^>=EI<^ s V>| 2 , 

s 

where Py r is the projection operator on the eigenmanifold V r corresponding to 
a r , and the sum S s is taken over a complete orthonormal set \(p^ > (s=l,2,...) 
of V r . The state of the system immediately after the observation is described by 
the vector Py r \i\) > . 

These assumptions are the main assumptions on which quantum theory is 
founded. In the canonical formalism they are used together with the assumption 
that the time evolution of a state vector is governed by the Schrodinger equation. 
In the formalism [3] this assumption is not used. Instead the assumptions (i) 
and (ii) are used together with the following postulate. 

(iii) The probability of an event is the absolute square of a complex 
number called the probability amplitude. The joint probability amplitude 
of a time-ordered sequence of events is product of the separate probability 
amplitudes of each of these events. The probability amplitude of an event which 
can happen in several different ways is a sum of the probability amplitudes for 
each of these ways. 

The statements of the assumption (iii) express the well-known law for the 
quantum-mechanical probabilities. Within the canonical formalism this law is 
derived as one of the consequences of the theory. However, in the Feynman 
approach this law is used as a basic postulate of the theory. 

It is also used the assumption that the time evolution of a quantum 
system is described by the evolution equation \^(t) >= U(t,to)\^(to) >, where 
U(t, to) is the unitary evolution operator 
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with the group property U(t,t')U(t',to) = U(t,to), U(to,to) = 1. Here we use 
the interaction picture. 

According to the assumption (iii), the probability amplitude of an event 
which can happen in several different ways is a sum of contributions from 
each alternative way. In particular, the amplitude < ^2|^(^o)|^i > can be 
represented as a sum of contributions from all alternative ways of realization 
of the corresponding evolution process. Dividing these alternatives in different 
classes, we can then analyze such a probability amplitude in different ways. For 
example, subprocesses with definite instants of the beginning and end of the 
interaction in the system can be considered as such alternatives. In this way 
the amplitude < ip2\U(t,to)\ipi > can be written in the form [3] 

< ^ 2 |£/(MoM >=< Hi>i > + l*dt 2 t dt x < ^|5(t 2 ,*i)Wi >, (2) 

Jto Jto 

where < ip2\S(t2,ti)\ipi > is the probability amplitude that if at time t\ the 
system was in the state \i[)\ >, then the interaction in the system will begin 
at time t\ and will end at time t 2 , and at this time the system will be in the 
state | -02 > • In the case of an isolated system the operator S(t2,h) can be 
represented in the form S(t2,ti) = exp(iHot2)T(t2 — t\)exp{— iHoti), Hq being 
the free Hamiltonian [3]. As has been shown in Ref.[3], for the evolution 
operator U(t,to) given by (2) to be unitary for any times to and t, the operator 
S(t2,ti) must satisfy the following equation: 

(*2-*i)S(t 2 ,*i) = tdt A PdtsiU-tdSfatJSfat!). (3) 

This equation allows one to obtain the operators Sfe, t\) for any t\ and t^, if the 
operators S(t f 2 , t[) corresponding to infinitesimal duration times r — t' 2 — t[ of 
interaction are known. It is natural to assume that most of the contribution to 
the evolution operator in the limit t 2 —> t\ comes from the processes associated 
with the fundamental interaction in the system under study. Denoting this 
contribution by Hi nt (t2,ti), we can write 

S(t 2 ,*i),-> H mt (t 2 M) + o(T e ). (4) 

Within the GQD the operator Hi nt {t^t-\) plays the role which the 
interaction Hamiltonian plays in the ordinary formulation of quantum theory: 
It generates dynamics of a system. Being a generalization of the interaction 
Hamiltonian, this operator is called the generalized interaction operator. The 
parameter e is determined by demanding that H int (t2, ti) must be so close to the 
solution of Eq. (3) in the limit £2 - * t\ that this equation has a unique solution 



If Hi n t(t2,ti) is specified, Eq.(3) allows one to find the operator S(t2,ti). 
Formula (2) can then be used to construct the evolution operator U(t,to) and 
accordingly the state vector \ip(t) >= \ip{to) > + // dh // 2 dtiS(t2,ti)\ifj(to) > at 
any time t. Thus Eq.(3) can be regarded as an equation of motion for states of 
a quantum system. By using (2), the evolution operator can be represented in 
the form [3] 

< n 2 \U(t, t )\ni >=< n 2 \ni > + 

_u — f°° A ex P[~ i( y z ~ E n 2 )t] < n 2 \T(z)\n 1 > exp[i(z - E ni )t ] 

+ 27ri-- (z - E n2 )(z - E ni ) ' [b) 

where z = x + iy, and y > 0, n stands for the entire set of discrete and continuous 
variables that characterize the system in full, \n > are the eigenvectors of the 
free Hamiltonian, and < ri2\T(z)\ni > is defined by 

roo ~ 

< ri2\T(z)\ni >= i / drexp(izr) < n2\T{r)\n\ > . (6) 

The generalized equation of motion (3) is equivalent to the following equation 
for the T matrix [3]: 



d < n2\T(z)\ni > < n2\T (z)\n >< n\T (z)\ni > . . 

dz = v ( z - E n y " () 

Thus, instead of solving Eq.(3), one can solve Eq.(7) for the T matrix. The 
operator T(r) and correspondingly the operator S(t2,ti) can then be obtained 
by using the Fourier transform 

~ l roo 

< n2\T(r)\ni >= — — - / dxexp(-izr) < n2\T(z)\n\ > . (8) 

At the same time, the T matrix can be directly used for constructing the 
evolution operator. According to (4) and (8), the operator T(z) has the 
following asymptotic behavior for \z\ — > oo : 

< ri2\T(z)\ni > — ► i / drexp(izr) < n2\H^(r)\ni > +o(|z|~^), (9) 

\z\ ^oo JO 

where f3 — 1 + e, and H\^ t {t2 — t\) = exp(—iHot2)Hi nt (t2,ti)exp(iHoti) is the 
interaction operator in the Schrodinger picture. As has been shown in Ref.[3], 
the dynamics governed by Eq.(3) is equivalent to the Hamiltonian dynamics in 
the case where the generalized interaction operator is of the form 

H mt {hM) = -2iS{t 2 -ti)fTj(ti), (10) 

Hi(ti) being the interaction Hamiltonian in the interaction picture. In this 



dl = —iHj(t)\i/j(t) > . The delta function 5(r) in (10) emphasizes the fact 
that in this case the fundamental interaction is instantaneous. At the same time, 
Eq.(3) permits the generalization to the case where the interaction generating 
the dynamics of a system is nonlocal in time [3,5]. 

As has been shown [3,5], there is one-to-one correspondence between 
nonlocality of interaction and the UV behavior of the matrix elements of the 
evolution operator as a function of momenta of particles: The interaction 
operator can be nonlocal in time only in the case where this behavior is "bad", 
i.e. in a local theory it results in UV divergences. In Ref. [6] it has been 
shown that after renormalization the dynamics of the three-dimensional theory 
of a neutral scalar field interacting through a (p 4 coupling is governed by the 
generalized dynamical equation (3) with a nonlocal-in-time interaction operator. 
This gives reason to expect that after renormalization the dynamics of an EFT 
is also governed by this equation with a nonlocal interaction operator. Let us 
now consider this problem on a toy model of the separable NN interaction. We 
will use the c.m.s., and denote the relative momentum of two nucleons by p 
and the reduced mass by \i. In this model the generalized interaction operator 
should be of the form 

< P2|#J2(T)|pi >= ^(P2MP1)/(T), (11) 

where f(r) is some function of duration time of interaction r. Let the form 
factor i/j(p) has the following UV behavior: 

Wp)~|p|~ a , (|p|-oo). (12) 

Assume, for example, that ip(p) = |p|~ a + #(?); an d in the limit |p|^oo the 
function g(p) satisfies the estimate g(p) = o(\p\~ 5 ), where S > a, S > |. In 
this case, the problem can be easily solved by using Eq.(7). Representing 
< p2|T(z)|pi > in the form < p 2 |T(z)|pi >= ip* (p2)ip(pi)t(z) , from (7), we get 
the equation 
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with the asymptotic condition 



t( Z ) r m+w), (14) 

\z\ — ► oo 



fi(z) = i Jq 00 drexp(izr)f(r), and Ek = #^. The solution of Eq.(13) with the 
initial condition t(a) = g a , where a G (— oo,0), is 



In the case a > h, the function t(z) tends to a constant as z — > — oo: t(z) — > A. 

z z^ — oo 

Thus in this case the function /i(z) must tend to A as z^— oo. From this it 
follows that only possible form of the function /(r) is /(r) = — 2zA(5(r) + /'(t), 
where the function f'{r) has no such a singularity at the point r = as the 
delta function. In this case the generalized interaction operator is of the form 
(10), and hence the dynamics generated by this operator is equivalent to the 
dynamics governed by the Schrodinger equation with the separable potential 
A^*(p2)^(pi). In this case the LS equation follows from Eq.(7), and our model 
is equivalent to the ordinary separable-potential model. 

Ordinary quantum mechanics does not permit the extension of the above 
model to the case a < |, because of the UV divergences. The GQD allows 
one to extend this model to the case — \ < a < \. In Ref. [5] such an extended 
model has been constructed for — \ < a < \. Let us now consider the case 
a — |. In this case the UV behavior of the form factor ^(p) gives rise to the 
logarithmic singularities in the Born series. At first we have to determine the 
class of the functions f\{z) and the value of j3 for which Eq.(13) has a unique 
solution having the asymptotic behavior (14). In the case a — |, the function 
t(z) given by (15) has the following behavior for \z\ — > oo : 

t(z) -> b 1 \n- 1 (-z) + b 2 \n- 2 (-z) + o(ln- 2 (-z)), (16) 

\z\ — ► oo 

where b\ = — (47r/i) _1 and 62 = &i ln(— a) — b\{M{a) + g^ 1 ) with 
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M(a) - / Irv jl ' ' d 3 k. (17) 

J a- E k 

The parameter 61 does not depend on g a . This means that all solutions of Eq.(13) 
have the same leading term in (16), and only the second term distinguishes the 
different solutions of this equation. Thus, in order to obtain a unique solution of 
Eq.(13), we must specify the first two terms in the asymptotic behavior of t(z) 
for \z\ ^00. From this it follows that the functions fi(z) must be of the form 

h(z) = h \n~\-z) + b 2 ln- 2 (-z), (18) 

and P — 2a — 1. By using (8) and the definition of fi{z), we can then obtain 
the functions f(r) and hence the generalized interaction operator 

< P2|^S( r )|Pi >= -^-^*(P2)^(Pi) l°° dxexp(-izr) I 1 + 2 2 

Ztt J -oo \m{—z) In \—z) 

(19) 

where z = x + iy, and y > 0. Thus, in the case of the UV behavior of the form 



the form (19), where for given ^(p) only the parameter 62 is free. The solution 
of Eq.(7) with this interaction operator is of the form 

< p 2 |T(z)|pi >= N(z)r (P2MP1), (20) 



where N(z) = b\ {-b 2 + 61 \n(-z) + M{z)b{] 

In Refs.[5] the model for — | < a < | (in this case the interaction operator 
is nonlocal in time) was used for describing the NN interaction at low energies. 
The motivation to use such a model for the parameterization of the NN forces is 
the fact that due to the quark and gluon degrees of freedom the NN interaction 
must be nonlocal in time. However, because of the separation of scales the 
system of hadrons should be regarded as a closed system. From this it follows 
that the dynamics of such a system is governed by Eq. (3) with nonlocal-in-time 
interaction operator. Note that the UV divergences of EFT's are the price for 
ignoring the quark and gluon degrees of freedom in describing hadron dynamics: 
EFT's are local theories, despite the existence of the external quark and gluon 
degrees of freedom. However, renormalization of EFT's gives rise to the fact 
that these theories become nonlocal. Below this will be shown by using the 
example of our toy model. 

As we have stated, the interaction operator (19) contains only one free 
parameter 62- However, if there is a bound state in the channel under study, 
then the parameter 62 is completely determined by demanding that the T matrix 
has the pole at the bound-state energy. For example, in the S S\ channel the T 
matrix has a pole at energy Eb = — 2.2246MeV. This means that \b(EB)]~ — 0, 
and, by putting a = Eb in Eq.(15), we get 



K*r-(*-*«)/* ViK-*) - (21) 

In this case 62 = &i hi(— Eb) — b\M(EB). Let us now show that renormalization 
of the LS equation with the separable potential leading to a logarithmic 
singularities produces the same T matrix. In [4] the problem of renormalization 
of the LS equation was considered by using the example of the separable potential 
< P2MP1 >= A^*(p2)V ; (Pi) wrtn VKp) = (d 2 + p 2 )~4. The corresponding T 
matrix is of the form < p2|T(z)|pi >= ^*(p2)V , (pi)^(^) with 

t(z) = (X- 1 - J(z))-\ (22) 

where the integral J(z) = / cftk z _^ has a logarithmic divergence. By using the 
dimensional regularization, one can get 



where J £ (z) = / d 3 e k ^_^ . The strength of the potential A £ is adjusted to give 
the correct bound-state energy: A^ 1 = J £ (E B ). In this way we get 

Uz)]- 1 = .L(E B ) - J,(z) = (z- E B )/rf^ __gWL_ (24) 

The right-hand side of (29) is well-behaved in the limit e — >• 0. It is easy to see 
that, taking this limit, we get the expression (22) for t(z). Thus renormalization 
of the LS equation with the above singular potential leads to the T matrix we 
have obtained by solving Eq.(7) with nonlocal- in-time interaction operator. This 
T matrix satisfies the generalized dynamical equation (7), but does not satisfy 
the LS equation. Correspondingly the Schrodinger equation is not valid in this 
case. The strength of the potential A £ tends to zero as e — > 0, and consequently 
the renormalized interaction Hamiltonian is also tend to zero. Thus despite 
the fact that the dynamics of each theory corresponding to the dimension 
T> = 3 — e is Hamiltonian dynamics for every e > 0, in the limiting case T> = 3 
we go beyond Hamiltonian dynamics. The dynamics of the renormalized theory 
is governed by the generalized equation of motion (3) with nonlocal-in-time 
interaction operator which in the case of the model under study is given by (19). 
To conclude, by using the model of the separable NN potential which 
gives rise to logarithmic singularities in the Born series, we have demonstrated 
that after renormalization the dynamics of a nucleon system is governed by the 
generalized dynamical equation with a nonlocal-in-time interaction operator. 
This gives reason to hope that the use of the GQD and a parameterization 
of the NN forces like (19) can open new possibilities for applying the EFT 
approach to the description of low-energy nucleon dynamics. By using an 
EFT, one can construct the generalized interaction operator consistent with the 
symmetries of QCD. This operator can then be used in Eq.(3) for describing 
nucleon dynamics. Note that in this case regularization and renormalization are 
needed only on the stage of determining the two-nucleon interaction operator. 
After this one does not face the problem of UV divergences. 
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